Permeability and conductivity of platelet-reinforced membranes and composites 
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We present large scale simulations of the diffusion constant D of a random composite consisting 
of aligned platelets with aspect ratio a/b >> 1 in a matrix (with diffusion constant Do) and find 
that D/Do = 1/(1 + cix + C2X 2 ), where x = avj/b and Vf is the platelet volume fraction. We 
demonstrate that for large aspect ratio platelets the pair term (x 2 ) dominates suggesting large 
property enhancements for these materials. However a small amount of face-to-face ordering of the 
platelets markedly degrades the efficiency of platelet reinforcement. 



Thin reinforcing platelets can be extremely effective at 
improving the barrier and in-plane mechanical proper- 
ties of composites and membranes. In particular there 
has been an explosion of interest in clay-reinforced ther- 
moplastics, thermosets and rubbers, with target applica- 
tions ranging from packaging to cars To achieve 
the theoretically promised enhancements requires well- 
aligned and well-dispersed clay platelets in these polymer 
matrices. 

The traditional theory of composite reinforcement is 
based on single inclusion theories which form a basis 
for self-consistent or effective-medium approximations. 
However the effect of aligned reinforcing platelets is not 
correctly described by single inclusion models [^[0], ex- 
cept at very low inclusion concentrations, even though 
many publications assume this approximation j^^] . The 
correct variable to use in describing platelet reinforce- 
ment in the large aspect ratio limit is the product of 
the aspect ratio (a/b) times the volume fraction (v/), 
x = avf/b. Since the aspect ratio of clay platelets ranges 
from 100 — 2000, x is typically not small as even a 1% 
inclusion volume fraction leads to large values of x. We 
calculate the diffusion constant in regimes where x is not 
small and derive a simple form which represents the data 
well. We find that the quadratic term dominates (i.e. 
a term proportional to x 2 ) and its dominance is due to 
narrow necks between platelets which are not included in 
single inclusion theories. 

However, chemically nano-dispersed clays tend to have 
face-to-face stacking, which leads to channels of matrix 
material through which diffusion is relatively easy. In 
a thermodynamic picture this corresponds to a phase 
separation of the material into platelet rich and platelet 
poor regions. We calculate the dependence of the dif- 
fusion constant on face-to-face alignment and show that 
the rule of mixtures works effectively. In clay-polymer 
materials deleterious platelet-poor channels are broken 
up by using extrusion or mechanical mixing. To date no- 
one has found a way to chemically modify the clays so 
that they prefer to arrange themselves in a more optimal 
(e.g. staggered) array. Nevertheless, theoretical studies 



suggest that nematic and staggered phases can be ther- 
modynamically stable JTof ] - 

The leading order term in the reduction in permeabil- 
ity due to platelet reinforcement is familiar in the com- 
posites community, though in a different context. There 
it is well known that a number density of aligned cracks, 
n = N/V (N is the number of cracks and V is the sample 
volume), reduces the conductivity of a material of initial 
conductivity (Tq according to, 
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for penny shaped cracks of radius a in a homogeneous 
three-dimensional medium fl_|| . The perpendicular (ie. 
normal to the crack surface) conductivity of cracked 
solids cr/ao, the perpendicular permeability, k/ko, of 
platelet-reinforced membranes and the perpendicular dif- 
fusion constant measured in platelet reinforced mem- 
branes, D/D , are related by, a/ao = k/ko = D(l — 
Vf)/D - However the experimental results for the per- 
meability of barrier films are presented as a function of 
inclusion volume fraction, v / . This is derived simply from 
Eq. (1) and (2) by using Vf — nv* , where v* is the vol- 
ume of an inclusion. Thus for barrier membranes, the 
leading order behavior for aligned rectangular sticks in 
two dimensions (v* — Aab) is, 
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and for aligned penny-shaped platelets of radius a and 
thickness 26 is, 
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The importance of the variable x = avj/b is evident from 
these expresssions. In the barrier film community the re- 
duction in permeability due to platelets is approximated 
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in a different way There it is argued that for high 

aspect ratio platelets, the increased tortuousity of typi- 
cal diffusion paths L p /L gives the qualitatively correct 
reduction in the diffusion constant, i.e. D/Dq ~ Lq/L p . 
The dependence of path tortuousity on inclusion volume 
fraction is given by L p /Lq = (1 + ai)f/2b), where a/b is 
the platelet aspect ratio. Note that this is of the same 
form as an effective medium theory based on equations 
(3) and (4). However, as we now show this tortuousity 
argument is qualitatively incorrect for platelet reinforced 
materials and if done correctly leads to a quadratic (ie. 
0((b/ (avf)) 2 )) reduction in the diffusivity in both two 
and three dimensions. 

We use the resistor respresentation to calculate the ef- 
fect of tortuous diffusion paths on the overall conduc- 
tivity, permeability and diffusivity. Consider a com- 
posite composed of randomly centered, aligned, non- 
overlapping, insulating sticks or pennies placed in a ma- 
trix of conductivity o~q. We define the typical perpen- 
dicular distance between inclusions to be I. The vol- 
ume fraction is related to I via, Vf w 26/(7 + 2b) « 2b /I 
as the inclusion volume fractions that are observed for 
high-aspect ratio materials are typically less than 10%. 
A tortuous path through a random array of these aligned 
platelets is approximated by a series combination of re- 
sistors each of which has typical resistance, 
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where po = 1/co- This resistance is calculated by consid- 
ering a "neck" of matrix material between two adjacent 
inclusions. This resistor has typical length of order a 
(we drop constant prefactors) and cross-section of order 
la d ~ 2 . In a film of thickness t, the resistance of a tortu- 
ous path is then R t = r t (t/(l + 2b)) « p ta 3 ~ d /l 2 . In a 
composite of transverse dimension L t (perpendicular to 
the thickness direction) however, there are many parallel 
paths of this sort and their conductances must be added 
to approximate the overall permeability or conductivity. 
The typical number of such paths is (L t /a) d_1 , so that 
the typical resistance of a composite of dimensions L d ~ 1 t 
is Rf w p ta 3 ~ d /((L t /a) d ~ 1 l 2 ). We are interested in the 
conductivity (cx permeability) which is related to the re- 
sistance by a = t/ Rf{L t ) d ^ 1 . Using Vf = b/l, we thus 
find 
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Note that the result (6) is due to the necks between 
platelets and is thus a pair term and is not included in 
single inclusion theories. In order to find typical val- 
ues for size of this quadratic term, and to compare its 
importance with the linear term found from single inclu- 
sion theories (Eq. (3) and (4)), we have carried out large 
scale simulations. To compare the above theory with the 
simulations, we need to include both the leading order 
term (Eq. (3) and (4)) and the quadratic term (6) in our 



analysis. The simplest form which contains both of these 
terms is, 
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where x = avf/b. As we show below this form works ex- 
tremely well for aligned non-overlapping, well-dispersed 
platelets over a broad range of concentrations. Impor- 
tant corrections occur if the platelets stick to each other 
or if the platelets overlap (ie. percolation effects) and we 
shall present the details of these effects elsewhere. 

Our calculation of the effective diffusivity of clay rein- 
forced membranes is carried out as follows. Clay platelets 
act as effective diffusion barriers to molecules such as 
oxygen and water, so we assign them diffusion constant 
zero. In contrast many polymers are oxygen and water 
permeable so we assign the matrix a finite diffusion con- 
stant Dq. We calculate the effective diffusion constant 
of a polymer containing volume fraction v/ of aligned 
platelets, as illustrated in Fig. 1. The effective diffu- 
sion constant of these composites is found by embedding 
the geometry in large square or cubic lattices. Random 
walkers are then started at random locations in these 
geometries. The end-to-end distance of these walkers is 
monitored as a function of time. Averaging the trajec- 
tories leads to the average diffusive behavior (r 2 ) cx Dt. 
The effective diffusion constant is then found by extract- 
ing the slope of a plot of (r 2 ) versus t. Our code for the 
"blind ant" jl3| method described above is very efficient, 
which enables simulations on very large lattices over long 
times. This is essential for the large aspect ratio compos- 
ites discussed here. Note that k/k — D(l — Vf)/D due 
to the fact that random walkers cannot be placed on the 
zero permeability barriers. In experiments, the perme- 
ability is measured by placing the barrier in a pressure 
gradient, with pressure p of gas (e.g. Oxygen) on one side 
of the membrane and the other side maintained at very 
low pressure. The steady-state flux of gas, /, through 
the membrane is measured and the permeability k = p/ f. 
The diffusion constant is measured by tracking the tracer 
diffusion of tagged particles, for example using NMR. 

A high precision test of Eq. (7) is presented in Fig. 
2a for aligned sticks in two dimensions and Fig. 2b for 
aligned squares in three dimensions. We plot the quantity 
(Dq/D—1)/x which, if Eq. (7) is correct, should be linear, 
ie. (D /D — l)/x = ci + c 2 x. The data of Fig. 2 confirms 
this form remarkably well over the entire range of inclu- 
sion concentrations, even close to the dense packing limit. 
A linear fit to the data of Fig. 2 yields the coefficients to 



the quadratic term, 
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0.165 ± 0.01 (from Fig. 2a), 



and c\ d = 0.050 ± 0.005(from Fig. 2b). The leading or- 
der coefficient in the concentration expansion, c±, is quite 
poorly converged in two dimensions. From Fig. 2a, it is 
seen that this c\ d is still decreasing even for large aspect 
ratio platelets. For the largest aspect ratios we studied 
its value is c 2d {a/b = 250) = 0.46 ±0.01. In three dimen- 
sions the convergence is much better and wc find c 3d — 
0.44±0.03. The theoretical values (fromEqs. (3) and (4)) 
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are ef * = tt/4 = 0.785..., and c^ enny = 4/(3tt) = 0.425... 
respectively. The numerical simulations are on lattices 
and so cannot be expected to be exactly the same as 
the continuum results (3) and (4). Nevertheless, it is 
clear that materials reinforced by well-dispersed platelets 
are not correctly described by effective medium theory 
based on a linear expansion in x, or the Nielsen formula 
(D/Dq — l/(l+x/2)) which is used in the diffusion com- 
munity ||[|. An important consequence of this result is 
that the property enhancements which are theoretically 
possible from platelet reinforced materials, for example 
clay-polymer nanocomposites, are much larger than has 
previously been suggested or observed. 

However the real morphology of, as synthesized, clay- 
polymer materials is illustrated in Fig. 3 The pres- 
ence of channels of matrix material (light regions in Fig. 
3a) and the strong face-to-face packing of the platelets 
is evident in Fig 3b. These materials do not provide 
the barrier performance promised by the result Eq. (7). 
The channels of matrix material act as a diffusion "short- 
circuit" so that the effective diffusion constant remains 
closer to that of the matrix material. In order to make 
the diffusion paths through the matrix material more tor- 
tuous, materials such as this are mixed in extruders or 
sheared in other ways in order to produce better platelet 
dispersion. The effect of various degrees of platelet phase 
separation on the diffusion constant is illustrated in Fig. 
4. The platelets are initially in a perfectly staggered ar- 
ray and then one sublattice of the array is shifted by an 
amount < s < a until the platelets are in a perfect 
face-to-face arrangement. The rule of mixtures works 
quite well for these phase separated systems (see Fig. 4), 
so that, 
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Clearly even a small amount of phase separation (finite s) 
leads to a significantly increased diffusion, that is the first 
term in Eq. (8) dominates for s > a/(c2X 2 ) — b 2 / (c^Vfa). 
In a very real sense, the larger the platelets are the more 
sensitive the system is to phase separation. For example 
if the reduction in diffusion constant promised by the for- 
mula (6) is 1000 fold, then any matrix channel, which is 
not highly tortuous, of size 1/1000 or larger will be dele- 
terious for the composite performance. This means that 
platelet dispersion must be very good indeed to achieve 
the full performance enhancements promised by large as- 
pect ratio platelets. 

In summary, we have shown that the effective diffusion 
constant (as well as the conductivity and permeability) 
of platelet reinforced composites is not well described by 
single inclusion theories, even though almost all exper- 
imental studies in the literature use expressions which 
are based on this limit. This is good news, as it means 
that the theoretically possible property enhancements are 
quadratic in the volume fraction, rather than linear. This 
general result should also apply to many other transport 



and mechanical properties of well-dispersed platelet rein- 
forced materials, and is due to the dominance of narrow 
necks in large aspect-ratio limit. However we also showed 
that face-to-face platelet ordering is extremely deleteri- 
ous to the performance of these materials (see Fig. 4 
and Eq. (8)), so that achieving the full enhancements 
promised by these materials remains a challenging syn- 
thesis problem. 
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FIG. 2. Tests of the form (7) of the text. A plot of 
(D /D — l)/x as a function of x = avf/b, where a/b is the 
platelet aspect ratio and Vf their volume fraction. The solid 
FIG. 1. Aligned square platelets randomly embedded in lines are fits to the largest aspect ratio data, (a) Data for 
a three-dimensional cube. The aspect ratio a/b = 25, the aligned slits randomly placed, without overlap, onto a square 
boundary conditions are periodic in all directions. lattice. The simulations were carried out on square lattices of 

size 2048 2 , over 6 x 10 4 steps in the blind ant algorithm (see 
text) and averaged over 30, 000 configurations, (b) Data for 
aligned squares (see Fig. 1) placed, without overlap, onto a 
cubic lattice. The simulations were carried out on cubic lat- 
tices of size 512 3 , over 5 x 10 4 steps in the blind ant algorithm 
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FIG. 3. Experimental platelet morphology at two length 
scales, (a) The upper figure shows that there is a phase 
separation into platelet rich and platelet poor regions, (b) 
The lower figure illustrates the face-to-face alignment of the 
platelets. The dark regions are clay, the light regions are the 
polymer (from ref. [14]). 
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FIG. 4. The effective diffusion constant as a function 
of face-to-face alignment of the platelets (of aspect ratio 
a/b = 128). The inset shows the geometry that was used for 
the calculations. Also included is the data for a random sys- 
tem with the same aspect ratio and it shows that the random 
system is similar to the perfectly staggered array of barriers, 
at least at the resolution of this plot. The three plots at inter- 
mediate s (dots, small dashes and large dashes) are fits using 
the rule of mixtures (Eq. (8)). 
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